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SUMMSiRY 


Numerical results of a theoretical investigation are 
presented to provide information about the effect of 
Variation of the diffei'ent design and opex-ating parameters 
on radial inflow turbine pei'formance . The effects of 
variations in tine mass flow rate, rotor tip Mach namJ:>ar, 
inlet flow angles, number of rotor blades and hub to 
shroud radius ratio, on tiie internal fluid dynamics of 
turbine rotors, was investigated. A procedure to estimate 
tlie flow deviation angles at the turbine exit is also 
presented aiid used to examine the influence of the 
operating conditions and tlie rotor geometrical configuration 
on these deviations. The significance of the results 
obtained are discussed with respect to improved turbine 
performance , 
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INTRODUCTION 


The benefits of using small radial inflow turbines in con- 
junction with automotive and aeronautical industries have been 
discussed by many authors [1, 2]. Boundary layer separation, 
secondary flow effects and heat transfer problems within the 
flow channels of these machines are examples of factors v^?hich 
impose a limit on their efficiency. These factors could br 
controlled by aerodynamic design methods based upon detailed 
knowledge of the flow behavior within the different passages 
of the turbine. Both theoretical and experimental investigations 
of flow conditions that result from various design configura- 
tions and operating parameters are considered necessax'y to 
develop this knowledge. 

In the present work a theoretical study, which makes use 
of recent advances in computational fluid mechanics, is 
presented. The object is to improve our understanding of the 
internal fluid dynamics of radial inflow turbine, rotors ujxder 
different operating conditions. A single stage rotor having 
a pressure ratio of 3,0 to 1 is used as a model in this stvidy. 
The effects of variations in the following design and operating 
parameters upon the rotor performance are investigated: 

a. Turbine flov/ rate 

b. Rotor tip Mach number 

c. Number of rotor blades 

d. Inlet flow angles to the rotor 

e. Hub to shroud radius ratio. 

The flow field analysis procedure and the computer program 
employed, to test the effect of the different parameters, are 
identical to that presented in Ref. [3] , Basically, the 
procedure is based upon the solution of the velocity gradient 
equation along an arbitrary quasi-orthogonal in the meri- 
dional plane. From this meridional flow solution the blade 
surface v'elocities are obtained according to the method 
presented in Ref. [4]. The factors which govern the selection 
of this procedure have been discussed in Pef. [5]. The 
assumptio.ns made in Ref. [3] yield to a solution of the flew 
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field which give good results only for the velocity distribution 
within the blade channels, away from turbine exit. 

Due to the peculiarity of the method of analysis of 
Ref. [3] in detecting the real flow behavior at turbine exit, 
an analytical model is developed in Part 2 of this study. 

The model offers a procedure whereby the flow deviation from 
blade surface at turbine exit is estimated. The model is 
based on the analytical solution of the governing aquatic^ ■ 
of two dimensional irrotational motion of flow between blades 
with some simplifying assumptions. The solution obtained 
allows also for the calculation of flow parameters at turbine 
exit in a quantitative manner, v/ith a more realistic dependence 
on the geometrical configuration and operating conditions. 

It is hoped that the results of the parametric study presented 
in Part 1, together with the application of the proposed model, 
yield an i.mproved design procedure for the radial turbine. 


PART I - FLOW ANALYSIS WITHIN BLADE PASSAGES 

General 

The detailed knowledge of the flow behavior through turbo- 
machines can be achieved through a field analysis of the 
unsteady, rotational, three dimensional viscous flow. One 
approach would be to attempt the solution of the complete Mavier 
Stokes equations. Clearly, such solution is out of reach nov; 
and in the immediate future. An alternative approach v'hich 
has been widely used is to assume that a two layer model is 
representative, i.e. an inviscid flov/ solution which interacts 
with an end wall boundary layer solution. 

Three dimensional inviscid flow analyses usually involve a 
combination of several two-dimensional solutions on intersactirg 
families of stream surfaces [5]. M.ost two-dimensional 
solutions are either on a blade-to-blade surface of revolution or 
on the mid-channel stream surface between tv?o blades (Fig. 1} . 

In t-he following study, the internal fluid dynamics of 
turbine rotors are investigated under different operating 
conditions. The viscous effects are neglected. The solution to 
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the equations of invlscid flow on the mid-channeX surface is 
obtained using the quasi-orthogonal method [3], This solution 
is termed here the "meridional flow solution". Using this 
solution, the blade surface velocities are obtained. 


Meridional Flow Solution 

The flov; field analysis procedure presented in Ref. [3], 
and employed here to study the flow behavior in the meridional 
plane of turbine rotors, is based on the streamline curvature 
method. In this method, the velocity gradient equation 
along a quasi-orthogonal is integrated numerically to calculate 
the distribution of flow properties. The governing equation 
for the velocity gradient is given as follows: 
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The coordinate system and the notations are shov/n in Figure 2. 

The numerical solution to equation (1) is assumed to be 
valid for the mean stream surface which extends from hub to 
shroud. This mean flow surface is considered to be parallel 
to the mean blade surface. Empirical corrections are made to 
take into account the difference between the actual flow angle 
and the blade angle at inlet to the rotor. Mo provisions 
have been reported in Ref. [3j to tai:e care of this difference 
at the turbine exit. In view of this fact, it is felt that 
it is quite meaningful to study the mechanism v/hich causes the 
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flew deviation from blade surface, at turbine exit, on a 
rigorous basis. In the second part of the present study, a 
calculation procedure will be developed to predict these flov/ 
deviations so that appropriate measures may be taken to 
improve the flow calculations. 


Blade to Blade Solution 

In order to render a rational judgement with regard to 
the effect of variation of the different design parameters 
upon the performance of the turbine rotor, knowledge of the 
velocity distribution over the rotor blades is necessax'y. 
With the velocities on the mean stream surface calculated 
according to equations (1) and (2) , the blade surface velo- 
cities can be obtained by several methods. One that gives 
good results, v/hen compared with the relaxation solution of 
the potential flov; equation, is given by Stantxz [4j. Ilis 
method is based on absolute irrotational flow and linear 
velocity distribution between blades. The follov/ing 
equations cased on these assumptions are used, herein, to 
obtain the required blade surface velocities. 


W 


COS3 coss ^ 

; — T ; — + u.' r ( t an S - t an B ) 

cos3 +cosb cosB s p 

s p p 


+ t(«r + Kslni) - tg)] 


(3) 


and 

W = 2W - W (4) 

P s 

In these equations, the velocity W and the flow angle B are 
those obtained from the meridional flow solution. The subscripts 
"s" and "p" refer to the suction and the pressure surfaces, 
respective ly . 
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C onditions of the Solution 

The flow behavior within the channels of a radial inflow 
turbine rotor is investigated under different operating conditions. 
The solutions v/hich will be presented are intended to show the 
effects upon the rotor performance of variations in the 
following parameters; 

1. Turbine flow rate 

2. Rotor tip Mach number 

3. Number of rotor blades 

4. Inlet flov/ angles to the rotor. 

The rotor considered in this study has a hub-shroud profile 
in the meridional plane similar to that shown in Figure 3a. 

Table 1 summarises the ranges over v/hich the different para- 
meters are varied. The standard solution presented in the Table 
corresponds to the following estimated operating conditions 
(design point of the turbine) : 


Turbine inlet total pressure, 

Turbine inlet total temperature, T ^2 

Turbine pressure ratio 

Rotational speed, r.p.m. 

Turbine mass flow rate, Q 

Stator nozzle exit flow angle, 

(inlet to the rotor) 

Working fluid 


= 42.5 Ib/in^ 

= 2710 “R 

= 3.0 

= 79 , 800 r.p.m. 

= 0.932 Ib/sec 

= 73.0 degrees 

= Air 


Other design point conditions needed for the analysis were 
determined theoretically from the previous data, and are given 
as follovAs: 


The angle between relative velocity vector 
and its meridional component at rotor 
inlet, 02 

Absolute fluid velocity at rotor inlet, 

Equivalent weight flow, (Q»^/6) 

Prewhirl at the rotor inlet, • r 2 

Additional drop in total pressure across 
the turbine rotor due to losses 


-40.91 degrees 
1575 ft/sec 
0.735 Ib/sec 
346.98 ft^/sec 
2.5 Ib/in^ 
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In order to obtain the reference conditions for the different 
operating cases of Table 1, a computer program is developed. 

The program requires as an input the stator and rotor geometries, 
the mass flow rate, the rotational speed, and the total 
conditions to the turbine. It calculates, through an iterative 
procedure using a one dimensional analysis, the absolute fluid 
velocity, V 2 , the angle 62 ^ and the prewhirl at the rotor 
inlet, V r« , These . alues are used as input data for the 

V 2 * 

computer program of Ref. [3] . A sample input to the program, 

for the turbine under consideration at the design point, is 

given in Table 2. The velocity diagrams at each mesh node 

(corresponding to the intersection of quasi-orthogonals with 

streamlines) j in the meridional plane, are also shown in 

Figures 3b and 3c for the design point conditions. 

For all cases studied and summarized in Table 1, the 

impeller passage is divided into eleven stream tubes. The 

inlet total pressure and temperature are kept constant at 

2 

values of 42,5 Ib/in and 2710“R, respectively. 

PRESENTATION OF RESULTS AND DISCUSSIONS 

Two groups of numerical results are considered in this 
section. The first one deals with the flow distribution in 
the meridional plane of the turbine rotor. The second group 
investigates the affect of variation of different parameters 
upon the velocity distribution over the rotor blades. 

Meridional (Hub-Shroud) Results 

The distribution of streamlines, lines of constant relative 
velocity and isobars contours in the meridional plane are 
considered. The solutions to be presented are selected to 
correspond no the cases 1, 4, 6 , 12, 14 and 3 of Table 1. 

a. Streamlines : The projection of the mean surface 

streamlines on the meridional plane for the different cases 
studied are shown in Figure 4. The streamlines are designated 
by a stream- function ratio such that the value of the streamline 
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indicates tht* percentage of the flow through the turbine between 
the strearolino and the rotor hub. The strearplino spacing is, 
therefore, indicative of the v'elocity relative to the rotor, 
with close spacing indicating high velocities and wide spacing 
indicating low velocities. It is obvious that the general 
trend of the streamline distribution is the sar.\e for all 
cases studied. 

b. Meridional velocity distribution ; Lines of constant 
relative velocity to the rotor are shown in Figure 5. For the 
design condition corresponding to the standard solution of 
Table 1, the rate of acceleration along the shroud is seen to 
be much higher than that along the hub (case no. 1). Such 
behavior is similar to that of flow approaching a turning duct. 
As expected, this rate of acceleration along the shroud is 
observed to increase with the increase of mass flow rate, 
rotational speed and the naiiber of impeller blades (compare 
cases 1, 6 and 12) , 

Examination of Figure 5 indicates that although the annulus 
cross-sectional area of the rotor increases by 73% from inlet 
to exit, the velocity increases generally all over the channel. 
This phenomenon .may be attributed ho the rapid decrease of 
density along the flow path as may be revealed from Figvure 6 . 

Along the hub, and especially near the rotor tip, regions 
of decelerated flow are observed for all cases studied. The 
e.xtent of these regions and the velocity gradient through 
each giv'e an indication of the boui’idary layer heha\'ior. This 
of course has serious effects or. the rotor efficiency. As a 
general rule, it is observed that for the cases of low mass 
flow the negative velocity gradient along the hub is so high 
that regions of reversed flow do exist, as shown in case no. 12. 
This would result in a high level of turbulence with 
accompanying flow separation and mixing losses. 

The present results have important significance on the 
process of channel design for a radial turbine, particularly 
when the rotor is operated at high specific speeds. In such 
cases, rhe i .nlet- to-exit radius ratro of the rotor is small, 
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moreover, the curvature of the hub is large to turn the flow 
from radial to axial direction. Thus a careful design of the 
rotor and flov/path is required to prevent flow separation 
from the hufe. 

Blade Surface Velocity Distribution 

The velocity distributions over the rotor blades are 
considered in what follov/s. The examples which are presented, 
illustrate the effect of changing the number of blades, the 
mass flov; rate and the impeller tip Mach number on these 
velocity distributions. This is follov/ed by an investigation 
of the effect of changing the stator exit angles upon the 
performance characteristics of the rotor. 

Effect of Changing Number of Blades : 

Plow analyses for the 10 and 14 bladed rotors, described 
in Table 1, v/ere also performed. Figures 7 and 8 compare the 
calculated blade surface velocity distributions for these 
rotors with the distribution of the basic 12 bladed rotor. 

At the design point (standard solution) the velocity distribution 
shouts some diffusion on the suction side of the blades 
especially at turbine exit near the shroud (Fig. 7) and at the 
middle channel near the hub (iig. 8) . The rotor has also a 
leading edge separation and flow reversal over the pressure 
surface near the hub and between the nondimensional distance 
from tip (Am/m) of 0.4 and 0.7 near the shroud. 

Examination of Figures 7 and 8 show clearly that an 
increase in the number of rotor blades is accompanied by a 
reduction in blade loading. This effect is expected since the 
increase in number of blades reduce the blade passage effective 
area through which the flow passes, giving more guidance to 
the flow. The last figures also indicate that reducing the 
number of blades result in an increase in the local diffusion 
as well as the extent of the flov/ reversal regions. 

It is a common practice to consider the pressure surface 
reversal regions as 'well as the suction surface diffusion as 
the governing parameters for rotor losses. Therefore, the 
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implication of the previous calculated velocity distribution 
is that a reduction of the number of rotor blades is 
accompanied by a loss of turbine efficiency. However, one 
should take into account that there ought to be an optimum 
for the total nmtber of blades. This optimum could be defined 
by weighing the opposing effect of increased friction losses 
with {or against) benefits gained by reducing the flow reversal 
region as a consequence of increasing the number of blades. 

Effect of Mass Flow Pate ; 

The effect of variation of mass flow rate upon the blade 
surface velocities is illustrated in Figures 9 and 10. In 
these figures, all design and operating conditions (other than 
the mass flow) were maintained constant at the standard 
value, (Table 1) . 

Figure 9 shows the calculated velocity distributions 
along the rotor hub at both the sucticn and the pressure 
surfaces. It is evident from this figure that the blade 
loading* increases as the mass flow passing through the turbine 
increases. Near the turbine exit section the increase in 
blade loading with mass flow results primarily from an increase 
in velocity along the suction surface. 

Figure 10 shows the calculated velocity distributions along 
the shroud of the rotor. At low mass flow the rotor tip element 
is operated ah an incidence angle above the design value and 
beyond the stall conditions as evidenced by the negative loading 
near the tip for the 80% mass flow case. The calculated 
negative velocity (in reversed flow regions) at tlie rotor 
pressure surface is noticed to be almost the same for all 
mass flow studied (reaching a maximum around 350 ft/sec) . On 
the other hand, the e.xtent of the region which the reversed 
flow occupies is inversely proportional to the mass flow rate. 

As for diffusion along the suction side, a relatively low 
value is obtained for the case corresponding to lov; mass flov.'. 


* Loading is defined here as the difference between the suction 
surface velocity and the pressure surface velocity. 
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The turbine rotor has only 510 ft/sec diffusion at a mass flow 
rate of BO‘4 of tlie design value comparod to 680 ft/sec for 
a 1004 mass flow rate. On the whole, it could be recognised 
from the last figures that the largest variation of impeller 
performance with mass flow occurred near the blade tip regions, 
thus indicating the criticality of design in these areas. 

This result has important impact on the process of blade 
design for new impellers , particularly when the design procedure 
is based on specifying the blade surface velocity distribution. 
The present results indicate that when high losses are 
expected near the inlet, the Loading due to blade turning should 
be small since substantial losses will occur due to rapidly 
decelerating flow near the rotor tip. 

Effect of Kotor Tip Mach Nuniber : 

The relative velocity distributions along the pressure and 
suction surface of the rotor blades, for different rotor tip 
Mach numbers are shown in Figures 11 and 12. In these figures 
all the design and operating conditions other than the rotational 
speed of the rotor were maintaii\ed constant (see Table 1). 
Basically, these figures indicate that a relatively high 
loading is initially formed near the rotor tip on both hub 
and shroud for low tip Mach numbers. On the other hand, as 
hip Mach Jn.ui\bers are increased, the loading near the impeller 
exit is observed to increase. The additional loading near 
the rotor tip may bo attributed to the fact that the rotor 
is operated at a relative inlet flow angle above it.s design 
value when the rotational speed is varied. These results 
demonstrate once more the sensitivity of the tLurbine rotor 
performance to the variation in inlet flow angles. 

.^s a general rule, it was observed at low r.p.ra. that the 
velocity distribution is severely distoi’ted. Not only the flow 
is reversed over a considerable portion of the pressure surfaces, 
but the deceleration along the suction surface is largo and 
will lead certainly to flow separation. The present results 
shows that high losses are expected near the turbine inlet, 
when the turbine operates below the design r.p.m. 
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Effect of Variable Stator Aiiqle on Turbine Performanc o: 

In the preceding discussions, it hes been demonstrated that 
the fluid dynamic characteristics of turbine rotors are highly 
sensitive to the variation in inlet flow angles. It is felt, 
therefore, that a detailed study of the effect of changing 
the stator exit angles upon the performance characteristics 
of the turbine will help to create a general uiiderstanding 
of the phenomena of rotor-stator interaction in radial 
machines . 

Three different stator configurations are considered here 
to investigate the effect of stator exit angles upon the mass 
flow characteristics of the turbine rotor. The meridional 
flow field analysis of Table 2 is repeatedly performed under 
different operating conditions of the rotor for each stator 
exit angle. Figure 13 shows the variation in equivalent flow, 
Q/W/6, with rotational speed and pressure ratio, 
across tiie rotor, for stator exit angles of 60°, 73° and 
80° respectively. In drawing these characteristics, the static 
pressure, p^, at exit from the rotor is evaluated on a mass- 
averaged basis. 

Figure 13 indicates that at pressure ratio above 3 and 
rotational speeds below 60% of design value, the turbine is 
choked with a limiting flow of 0.827 Ib/sec for the 73° stator 
exit angle. At all other conditions, for the same stator exit 
angle, weight flow decreases rapidly with increasing speed, 
indicating the impact of centrifugal effect upon the turbine 
rotor performance. It may be noted also that slight changes 
in pressure ratio have a large effect on weight flow at all 
rotational speeds. 

Comparison of the three sets of curves of Figure 13, 
which correspond to different stator exit angles, shows clearly 
the large influence of stator exit angles upon the choking 
conditions in the radial machine. As a matter of fact, the 
choking mass flow is reduced by about a third v;hen the stator 
exit angles are reduced from 80° to 63° at a rotational speed 
of 100% of the design value. 
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PART II 


FLOW DISTRIBUTION AT TURBINE ROTOR EXIT 


PRELIMINARY CONSIDERA.TIQMS 

Experience has shown that the shape and the geometrical 
configuration of the turbine exit section has a serious effect on 
its performance [6]. It has been common practice to design 
turbine rotors such that the flow leaves the machine in an 
axial direction. The exit section of the rotor is often 
designed under the assumption that the relative flow leaves 
the blade tangent to its trailing edge. In fact, the relative 
exit flow angle, 6^, at each radius of the turbine exit, 
usually deviates from the corresponding geometrical exit angle, 

3 , (Fig. 14). Such deviation may be attributed primarily 
to the secondary flow effects caused by the relative swirl 
motion within the various blade passages due to the impeller 
rotation, as shown in Figure 15. An approach v/hich has been 
used to determine these deviation angles is based upon the use 
of an empirical formula [7]. As nev/ design concepts are 
developed, a need for a more comprehensive and rigorous method 
of predicting the deviation angles is needed. This step is 
critical in the design process since excessive flow angle 
deviations may result in flow separation with the subsequent 
drop in the efficiency of the expansion process. 

In the following study, an analytical model which predicts 
the magnitude of swirl velocity behind the rotor of a radial 
inflow turbine is proposed. The model is based on the analytical 
solution of the governing equations of two dimensional 
irrotational flow between blades with some simplifying 
assumptions. The effects of variations in the following design 
parameters of the impeller upon the swirl velocity distribution 
are also investigated: 

a. Hub to sliroud radius ratio 

b. Circumferential spacing between blades 

c. Impeller rotational speed 

d. Turbine flow rate. 

From the analytical solution obtained, along with appropriate 
assumptions concerning the blade shapes, deviation angles ab 
rotor exit are calculated. In this context, the flow angle 
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distribution in tangential direction, along the turbine exit 
section is presented for the rotor of Figure 3a under various 
operating conditions. Finally, possible improvements in flow 
non uniformity by means of design changes using the results of 
the present method analysis are indicated. 


ANALYSIS 

Method of Approach 

The method of approach taken here is to construct a 
calculational procedure for analyzing the flow conditions at 
the rotor exit. The analytical model developed for that reason 
is based upon the analysis of two dimensional flow between 
blades on the flow surface of revolution, S , of Figure 16. 

A developed view of the flov/ stream channel, of thickness, h, 

in the meridional plane, together with the coordinate system 
used in the analysis are also shown in the same figure. The 
flow field aiialysis procedure employed to derive the necessary 
governing equations used in the present study is identical to 
that presented in Ref. [8]. The variations in the flow 
properties between blades on the flow surface of revolution, 
is obtained through the solution of the governing 
equations by the method of sep.aration of variables. The 
solution obtained will be used to predict the influence of 
the chosen blade shape {or boundary conditions) upon the 
relative velocity and the flow angle distribution at exit 
from the turbine rotor. > 

Governing Equations 

a. Continuitv: 

— 

From continuity consideratidi”is for a steady compressible 
flow through th.j stream tube of Figure 16, the following 
expression is obtained: 

3^ (j_ V H R) + U H) = 0 (5) 

o o 
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A nondimensional stream function, tp, satisfies the continuity 
Equation (5) if defined as 



V H P. 


(5a) 



v/here 
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V 




H = 



M = 


TO 
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(5b) 


(5c) 


b . Relative Irrotational Motion 

In the absence of shocks and viscous effects, the absolute 

motion of the fluid on the surface of revolution, S , is 

m 

irrotational. Therefore the relative circulation is zero. 
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Substitution of the stream function, i)t, as defined by Equation 
(5a) and (5b) in the above equation gives 
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where 


sina 


dR 

dM 


Equation (7) is a nonlinear partial differential equation for the 
nondimensional stream-function distribution of two-dimensional 
flow past arbitrary blade shapes. For given boundary conditions, 
equation (7) could be solved using numerical techniques to obtain 
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the velocity field between blades [9]. Basically, these numerical 
solutions require a large amount of computation time. Therefore, 
it would be advantageous to have qvticher means, based on some 
simplifying assumptions, for estimating the flow conditions 
between blades. With this motivation, a mathematical solution 
to the governing equation is considered in the next section. 


Assumptions and Limitations for the Analytical So3ution 

To simplify the task of finding an analytical solution to 
equation (7), the following assumptions have been made: 

a. The flow through the passage betv/een blades is 
considered to be incompressible. 

b. The passage height in the meridional plane is 
constant with radius, hence H = 1. 

c. The rotating impeller has straight radial blades. 

d. Variation of axisymme trie stream surface angle, ct , 
with radius is neglected. Consequently, the flow 
surface of revolution is conical, and a is constant. 

Under the stipulation of the foregoing, equation (7) is simplified 
to give 


2 


sin a 



sin a 
P 3R 


9 

1 

2 . 2 
R sinci 30 


(8a) 


also, equations (5a) and (5b) reduce to 


V 


1 

R 30 


U 


? 1*1 
3R 


(8h) 


Equation (Sa) is a Poisson equation and its solution is given by 
the method of separation of variables as; 

R(P.) (9) 

where ' in the above equatica satisfies Laplace equation, then 


T|; = (R) 3(0) 4- (P) 


( 10 ) 
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Substituting for f in equation {8) one obtains s 

dR, 


,2 n d^R, dR, , .2- 

sin tt 1 , _ 1 1 , 1 3 9 _ rt 

-rT” ^ dT^ ^ 0 :t2 " ° 


dR 


39 


UOa) 


Similarly, substituting for R 2 (R) in equation (8) one obtains: 


*^"^2 1 ^ „ _2 ^ Q 

^j^2 R R sinct 


UOb) 


Therefore, the expressions given by equation (10) satisfies the 
governing equation (8) providing the following relations are 
satisfied; 


R — ~ + R 
dR^ 


dR^ 

dl" 


R. 


+ 


. 2 
sxn a 


0 


1 d" 9 

e d9“ 



d^R^ 

dR“ 



_2 

sina 


0 


(11) 


where \ = constant, to be determined from boundary conditions. 

After solving the set of the above given ordinary differential 
equations, the general solution of the stream function is 
as follows; 






+ C, 


in R) + A-sin(— r — in R) ] [C, sinhxe 

M «!|} -U 11 ^ X 


( 12 ) 


where A^, A 2 , C^, C., , and are constants of integration to 

be determined from boundary conditions. 

Once the value of is deduced from equation (12) , the 
complete analysis of the theoretical velocity distribution in 
the rotating impeller channel is obtained using equation (8b) . 
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If the angular spacing between two successive blades is denoted 

by 2a, and the radius ratio of the hub (= ^hub'^^shroud^ ^1’ 
The boundary conditions to equation (12) in the special case 
of a pure radial turbine, for which the bl^xdie. Jslexole si'reexw 
surface is similar to that given in Figure 16C, may be written 
as : 


1. 

O 

II 

at 

0 = 

± a 

2 . 

II 

o 

at 

R = 

1 (corresponding to shroud radius) 

3. 

if; = 0 

at 

R = 

^1 

4 . 

•Me) = -Tj; (0) 

at 

any 

radius ratio. 


The first boundary condition gives; 


The second boundary condition results in; 



^2 ~ 2 s in a 

VJhile the third boundary 

^ _ n-r sina 
' ” ^'l 

- 1 ) 

D = - — ^ 

1 2 in rj^sina 

Hence, the expressio 
reduces to 


(For Laplace equation) 

(For Poisson equation) 

condition gives: 

(For Laplace equation) 

(For Poisson equation) 

■n for using the superposition 


principle , 


V * / inP» i,n’TSinciY^’i 

) sin{nT ) cosh (— — )e 

n in r in r- 

n=l 1 1 


2 sine 




(13) 
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In order to determine the coefficient (A^) in equation (13) 
the first boundary condition is applied once more: 


I sin(n,r 1^) costl ( 


n=l 


nn sxncti 

JKF. 


2sina 




( 14 ) 


Expanding the right hand side of equation (14) in a Fourier 
sine series, one obtains: 




niT 


^n = 


nr +4 (An r^) 


[rsin«] [cosh ] 


(15) 


The general solution for may now be v/ritten as 

^ -1 

1 n=l 



mr 


2 2 
n TT 


+ 4 ( in 


-] [- 


sin ( 


nir in P , 
in ri 


ri> 


IT sina 


) cos 
-] [ 




sing 

^1 


) e 




(16) 


It is to be noted that this latter solution is obtained from 
the irrotationality condition. Then, che stream function i|j 
expresses the case of a displacement flow, that results from 
the rotation of blades in still fluid. The pure through-flow 
v;hich accounts for the net discharge of the impelleir is obtained 
in a traditional way by placing a line sink at the center of 
rotation of stationary blades. The nondimens ional stream function 
for a through-flow at a volume flow rate Q may be written as 
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(9 + o) 


(17) 


’f' 5 

2n ojrj h_ 
s s 

Various operating conditions may be obtained by a linear com- 
bination of solutions (16) and (17). This procedure is valid 
by the linearity of Laplace equation and the fact that each 
solution separately satisfies its own boundary condition. 


ESTIMATION OF THE PELATIVE SWIRL VELOCITY DISTPIBUTIOM 

Shown in Fig, 17 is the blade exit configuration in a typical 

turbine. The geometrical blade angle (0 ) is related to the 

5 

minimum width of the flow passage at blade exit (o) and the blade 
pitch (t) at certain radius by the relation 

g = sin ^ (o/t) (18) 

In the absence of secondary flow caused by sv;irl motion, the 
relative flow at the blade exit would have a direction of 0 . 

When the effect of relative swirl velocity, u , is considered, 
the relative exit flow angle, 0^, at each radius of the turbine 
deviates from the corresponding geometrical angle, At the 

shroud, swirl motion exhibits backward characteristics because 
the fluid actually leaves with a slight backv;ard component, 
as illustrated in Fig. 18a. At the hub, the relative swirl 
motion produces a forward characteristics as shown in Fig. 18b. 

The important effect, from the practical point of view, is that 
the tangential component of the absolute velocity, is 

altered across the rotor exit in accordance with the sv/irl 
velocity distribution. Consequently, radially nonuniform exit 
flov; conditions prevail with the subsequent redistribution of 
available power output at each radius of the turbine exit. 

In order to obtain an accurate estimation of the net power 
output from a turbine for design purposes, a knovjledge of the 
relative sv;irl velocity dis ti'ibution at rotor exit is required. 
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The task of producing the necessary equation suitable for the 
prediction of relative swirl velocity is presented below. 

Referring to Figure 18, the tangential component of swirl 
velocity is given by; 

u = cor - 


(19) 


From the expression given by equations (8b), (16) and (19), 

the value of the nondimensional tangential component of swirl 

velocity, U , at any radius, R, will be; 
s 


= R - sin a 
s 


DR 




. 1,^1 ‘ 
'iRHn r 




n=l 




TT Rsina 


,X2.nR^ 

cos (rn;~) 


cosh (\0 ) 


sin a cosh (Xc ) 


(20) 


v/here 

u 

s tor 

s 


_ n*r since _ jr 

in r^ ~ z 


r 


[1~ (1)^^ r^^] 


nir 

■J 

n“Tr“ + 


9 J 


4 (in 1 - ) 


(2Qa) 


The theoretical distribution of (U^) as given by aquation 
(20) is seen to depend primarily upon the follov;ing parameters: 

1. Number of blades. 

2. hub to shroud radius ratio, r^. 

3. The radius (R) and position (6) at which swirl velocity 
is evaluated. 

4. Axisymmetric stream surface angle, a, independent of 
the turbine flow rate. 

Variation of v/ith these parameters will be studied in the 
next section, for the special case of a radial inflow' turbine 
in w’hich, a = ;i/2 . 
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IV. RESULTS AND DISCUSSION 


TS'lo groups of numerical examples are v/orked out using the 
present method of analysis. The first one investigates the 
effect of number of blades and hi.ib- to- shroud radius ratio on 
the swirl velocity distribution at turbine exit. The second 
group deals with the influence of the same parameters in 
addition to the effect of mass flow rate and the impeller 
rotational speed on the flov; angle distribution at rotor 
exit. 

Passage"Averaged Value of Swirl Velocity Distribution at 
Turbine Exit 

Calculations which give an indication of the magnitude of 

the tangential component of swirl velocity, U , at turbine exit 

s 

are performed using the first ten terms of the series in 

equation (20). Higher order terms have been neglected. Sample 

results showed that the maximum error introduced by such 

approximations lie within 0,2% as compared v/ith those calculated 

using up to fifty terms. The value of U from blade-to-blade 

s 

at different angular positions, 0, are averaged over the blade 

passage at each radius ratio, R. The variation of this average 

value, U , near the hub is plotted in Figure 19 as function of 

hub to shroud radius ratio, r^, with the numl^er of impeller 

blades, Z ^ as a parameter. Absolute value of the swirl velocity, 

U^, is seen to decrease considerably as the number of blades 

is increased for all values of r^ . These curves also indicate 

that for a given nuniber of blades, Z, increasing the hub to 

shroud radius ratio, increases the absolute value of U until 

s 

a maximum is reached. The points of maxima are shifted towards 

lower values of as the number of blades decreases. 

Figure 20 shows the variation of U near the shroud. 

s 

Unlike the previous case, it is obvious that for a given number 
of blades, there exists a region of relatively constant U , 

« * 3 

and a region of changing . The region or constancy extends 
over a large portion of radius ratio r^, and its extent increases 
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as the nuiT±)er tf blades increases. Accordingly, a reduction 
in the absolute value of U at the shroud is achieved only by 
selecting higher design values of hub-to-shroud radius ratio. 

The effect of variation of both the hub- to~shroud radius 
ratio, and the angular spacing between blades, a, upon 
the swirl velocity distribution at turbine exit is illustrated 
in Figures 21 through 23. 

Figure 21 shows the radial distribution of the mean value 

of U over the exit section of the bladed impeller having 
s 

different design values of r^. It is evident from this 
figure that the hub-to-shroud radius ratio has a serious 
effect on the U distribution. For lov/ values of r, , the 
distribution obtained is nonlinear from hub to mid radius 
approximately. However, the distribution becomes generally 
linear from hub to shroud as the radius ratio approaches 0.7. 
Figure 21 also shov;s that the center of rotation of the swirl 
motion in an impeller having equal to 0.2 occurs at R = 0.7 35, 
and for a 0.7 radius ratio impeller at 0.845. This indicates 
that the radius of rotation of fluid decreases as the hub 
to shroud radius ratio of the impeller decreases. 

The effect of changing the angular blade spacing, 2a, 

{by varying the nurober of blades) on the radial distribution 
of is shown for the case of a 6 and a 10 bladed impeller 
in Figures 22 and 23 respectively. The variation of U with R 
exhibits the same trends as discussed for Figure 21. Hov?ever, 
the magnitude of swirl velocity, , is decreased by increasing 
the number of blades. 

The impact of increasing the number of blades, upon the 
position of center of rotation of fluid at turbine exit 
is also shown in Figure 23. It is evident from these 
figures that the center of rotation shifts towards the shroud 
as the number of blades is increased. 
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Flow Angle Distribut3.on at Rotor Exit 

In studying the effect of swirl velocity, U , upon the 
relative exit flow angle, 3^^ from the turbine rotor, both the 
flow coefficient, cji , amd the geometrical angle distribution, 

3 , at turbine exit have to be known. The flow coefficient 
is related to the meridional component of absolute velocity, 
t of Figure 18 by 

= ^3 

*^^shroud 

and is directly proportional to the flow rate passing through 
the rotor. Representative values of used in the present 
study are chosen to correspond to the Cases 1, 2, 3. 7, 8 
given in Table 1, The relation between impeller exit radius 
and geometric exit angles used to carry out the necessary 
calculations is shown in Figure 24. 

Sample results showing the influence of swirl velocity 
upon the distribution of relative exit flow angles are shown 
in Figures 25 and 26. These figures indicate clearly that 
the flow loses the directional influence of the blade surface 
for a finite number of blades. 

E’or a given impeller geometry, hhe difference between and 6 
depends upon the relative magnitude of the mass flov/ rate and the 
impeller rotational speed- These results indicate that for 
high specific speed, rotors, as in the case of small radial 
turbine, large deviation angles can occur even for sizable 
values of the flow rate. 


* Deviation angles are defined here as the difference betv;een 
relative exit flow angle, 33 , aiid the corresponding geometrical 
angle (A = 
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For all cases studied, a positive deviation of flow angles 
over the geometrical angles are observed near hrib, while the 
flow near the tip tends to have locally negative deviation 
angles. This implies that near the hub, separation of flow 
is liable to occur at blade suction surface, on the other hand 
near the impeller tip separation may take place at blade 
pressure surface. This conclusion conforms with the experi- 
mental work presented in reference [10]. 

The effect of increasing the number of blades is to 
decrease the predicted deviation at all radii, Fig. 26. The 
position of maximum deviation is seen to occur near hub for 
low mass flow, while it shifts tovsrards the shroud for high 
mass flow rates. Fig. 25. 

Briefly, it may be concluded that the effect of swirl 

velocity, U , is to produce a non-symmetrical devicition angles 
s 

distribution at rotor exit. The deviation angles at all radii 
may be reduced by increasing the flow rate, decreasing the 
impeller rotational speed, decreasing the angular blade spacing 
and choosing an optimum hub-to-shroud radius ratio. 

It is obvious that in order to overcome the adverse 
effects of non-uniform flow conditions at tur-bine exit, 
procedures with which the effects of appropriate counter- 
measures can be estimated must be available. It is felt that 
with the present theoretical model such estimation can be 
adequately accomplished. 


25 



REFERENCES 


1. Wood, Homer J. , "Current Technology of Radial Inflow Turbines 
for Compressible Fluids," Trans. ASME, January 1963, p. 12. 

2. Wallace, P.J,, "Theoretical Assessment of the Performance 
Characteristics of Inward Radial Plow Turbines," Proc. 

Inst. Mech. Engrs . , 1958, pp. 172, 931. 

3. Katsanis, T. , "Use of Arbitrary Quasi-Orthogonals for 
Calculating Plow Distribution in the Meridional Plane of a 
Turbomachine," NASA TN D-2546, 1964. 

4. Stantiz, J.D. and Prian, V.D., "A Papid Approximation Method 
for Determining Velocity Distribution on Impeller of Centri- 
fugal Compressors," NACA TN 2421, 1951. 

5. Wu, ChLuig-Hua, "A General Theory of Three Dimensional Flow 
in Turboraachines of Axial, Radial, and Mixed Flow Types," 

NACA TN 2604, January 1952. 

6. UACL Engineering Report No. 458, "Phase III Interim Report 
No. '6, 90 Degree liiward Radial Turbine Research Program," 

Project P3. 

7. Mizuraachi, T. and Kitano, M. , "Study of Aerodynamic Character- 
istics of Rotating Blades in a Radial Inflow Turbine," paper 
presented at Tokyo International Conference, Oct. 4-8, 1971. 

8. Vavra, M. , "Aero thermodynamics and Fluid Flow in Turboniachinery , 
John Wiley, Nev; York, 1960. 

9. Katsanis, T. , "Computer Program for Calculating Velocities 
and Streamlines on a Blade-to-Blade Stream Surface of a 
Turbomachine," N.ASA TND-4525 , 1968. 

10. Lakshminarayana , B,, "Three Dimensional Flow Field in Rocket 
Pump Inducers," Transactions of the ASME, Journal of Fluids 
Engineering, December 1973, pp. 567-577. 


26 



NOMEMC’l.ATUPi: 


Constants of integration in equations (11) and (14) . 
Absolute fluid velocity, ft/sec. 

Stream channel thickness normal to meridional stream 
surface, ft., or static enthalpy, ft. lb. /slug 

Normalized stream-channel thickness defined as h/hg 

Distance along meridional stream surface, ft. 

Normalised meridional stream line distance defined 
as m/r 

s 

Throat dimension, ft. 

2 

Total pressure, Ib/in 

Mass flow passing through the turbine, Ib/sec. 

Radius from axis of rotation, ft. 

Normalized radius ratio, defined as r/r 

s 

Distance along arbitrary quasi-or thogonal in meridional 
plane 

Blade thickness in tangential direction 

Tangential component of swirl velocity, ft/sec. 

Normalized tangential component of fluid relative 
velocity, defined as W^/wr^ 

Normalized tangential component of swirl velocity, 

defined as u /wr„ 

s 

Normalized meridional component of fluid relative 
velocity, defined as W /wr 

Fluid velocity relative to blade, ft/sec. 

Prerotation in equation (1) , ft*^/sec. 


Axial coordinate 



z 


Number of blades 


ct Aiigle between meridional stream line and axis of 

rotation, deg. 

6 Angle between relative velocity vector and meridional 

plane, radians 

r Circulation 

«S Ratio of turbine inlet total pressure to atmosphere pressu 

9 Relative angular coordinate, deg. 

\ Coefficient in equation (11) 

p Mass density, slugs/cu.ft. 

Flow coefficient at rotor exit, defined as C, /er 

Im s 

^ Normalised stream function 

3 

p Density, lb/ ft 

2a Angular blade spacing defined as (2rr/Z) 

ui Rotational speed of the impeller, rad/sec. 

Subscripts 

2 Inlet to rotor 

3 Exit from rotor 

o Total conditions 

g Geometrical 

s Shroud 

9,r,ra,x Components of argument in tangential, radial, 

meridional, and axial directions, respectively. 

Superscripts 

Mecin value 
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TABLE 1 


PARAl'lETERS FOP NUMERICAL SOLUTIONS 


lution 

Mo. of 
blades 

Nozzle angle 
(degree) 

Mass flov; rate 
(% of design) 

Rotational speed 
{% of design) 

Rotor inlet 
relative sv;j.rl 
angle, 

Prev;hirl 
at rotor 
inlet, 
ft^/sec 

1 

12 

73 

100% 

100% 

-40.91 

346.98 

2 

10 

73 

100% 

100% 

-40.01 

346.98 

3 

14 

73 

100% 

100% 

-40.91 

346.98 

4 

12 

73 

100% 

120% 

-61.93 

346.98 

5 

12 

73 

100% 

80% 

- 2.2361 

346.98 

6 

12 

73 

100% 

60% 

38.25 

346.98 

7 

12 

73 

90% 

100% 

-59.61 

288.49 

8 

12 

73 

80% 

100% 

-69.47 

241.59 

9 

12 

60 

160% 

100% 

-33.54 

300.62 

10 

12 

60 

160% 

80% 

-16.22 

300.62 

11 

12 

60 

160% 

60% 

12.12 

300.62 

12 

12 

80 

65% 

100% 

23.60 

485.61 

13 

12 

80 

65% 

80% 

57-51 

485.61 

14 

12 

80 

65% 

60% 

68.93 

485.61 


Remarks 


Standard 



2; INPUT C?.PD LISTING ORlGlNAi; PAGE IS 

OF POOR QUALITYI 

a) Inlet Flow Conditions and Rotor Geometry Data; 


10 11 17 13 ft3b0;6fen? Q.9-3g.a m.OOOO 1.33 00 1715.0000 

0000 3710,0000 3«b.9Ro5 ' 0,0«20 0.0010 2.5000 

_1_3 1 3 1 1 . OO - OO «4-r^WK>-0 0 . 9 1 1 .8^ 00 

0,00001 

0.5000 0,5250 0,5750 0,6500 0.7250 0.R125 0.8B75 

l-r0*^>-3-§ l-r-W>2-5 — — -- 

0,0 0,0093 0.0250 0.0750 0,1625 0,2625 0,0125 

0.5875 0.8750 1.1250 

_2_,J7jXOO 2-^V^ 2-T-5WI0 2V-3-6-5-5 2.3500 2—1^7^^ 

2.0750 2,0500 2.0375 

2,7000 2,5000 2,2500 1.0750 1,5500 1,3000 1.0250 

— Q—-g-l-g-5 0 , 53 7 5 0 ^ . 5 0 0 0 — — "• ~ 

0.0 0.0 0,0 - 0 , 000 a -0,0027 -0,0090 -0.0200 

-0.0517 -0.0972 -0,1632 -0,2087 -0,3512 -0.0660 


b) Blade Thickness Table; 


-G-jrO ^ — 0 -r- t-Q 0 n 0 -- 3 . [} ^ 0,3000 0 - ; 0 000 0-r5“0-e-0 O. b OOO - 

0.7000 0,8000 0,9000 1.0000 1.1000 1.2000 

0,3850 0.3580 0.3220 0.2900 0.2600 0,2300 0.2010 

-0— 1^^^-e O-r-l-5-e-O O-_3-04Vf) 0 . Q9g0 

0,3750 0.3050 0,3110 0.2800 0.2500 0.2200 0,1910 

0.1 650 O.IOOO 0 , 1 1 90 0 , 1 000 O.O8OO 0,07iJO 

-G-r-;W>5-0 <H 3 3 3 0 0-;-34)-0O O-r^^O-O (>^-2-0-!^! 0-r2-HH) 0-^HM-O- 

0,1550 0.1310 0.1100 0.0920 0.0780 0.0690 

0.3550 0.3210 0.2900 0.2590 0.2290 0,2000 0.1710 

4)-,-l^ix0 O^-l-2-l-0 G-rO-7-3-^) 

0,3050 0,3100 0.2800 0,209 0 0,2190 0.1890 0.1610 

0,1360 0,1130 0.0930 0,0790 0.0690 0.0590 

^-.Ji_3-3^3 0-,-^-0-OJ3 iV2_f>O0 0-^2-3-a.O y^JsiQ 

0.12?0 O.iOilO 0,0860 0.0730 0.0 630 0,0560 

0,3200 0.2900 0.2580 0.2270 0.1970 0.1680 O.UOO 

■0^ 11 7 0 5-9 u_a^a9 (VO-6 h?.j> 9-yj;-5-Q-r 

0.3100 0,2790 0.2470 0.2150. 0.18^0 0.1570 0.1300 

0.1070 0.0870 0.0730 0.0620 0.0550 0.0090 

4)-._34MW2 Ot- 2-^^8-0 ^ Ch,-2-97UI — « 0-r4-‘ULO 0. 12Q -&- 

0,0980 0,0790 O.O680 0.0580 0.0510 0,0060 
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Blade Thickness Table CContinuad) 


0,2«00 

0.2510 

0.0790 

0,0700 
0, 1940 
-o^-O-b-txa- 
0,1 660 
0.0500 

0,0^00 
0 , 1 n 9 0 

^^-Q3-0.0_ 

o,o«oo 

0,0200 

■0 . 0 

0.7000 
0 . 5 0 0 n 
.l-.iiAao- 
2,4600 


0,2570 

_n-j 4 l 3 Lj.O 

0 , 2^160 

0,0670 

O.OfalO 
0. 1940 
O'S 0 0 ~ 

0 , 1 660 
0,0900 
_O.,-U-7-0_ 
0 , 0300 
0. 1090 
_a-,-o? 4 ).o — 
0 , ofloo 
0,0100 
_fl.,.u 3 .ao — 
0 . flOOO 
0 ,6400 
— V -,-6 2 - 0 - 0 — '■ 
2.6000 


0.2240 

0.2130 
0,0570 
-tt-r 24 V?.(V- 
0,0520 
0, 1900 

0.1660 

0.0300 

-O^UVLO- 

0.0200 

0,1090 

-O-.-OlJUCL- 

0,0600 

U.O 


0.1930 
- 0 -J- 0 J 544 )- 
0,1820 
0,0500 
- 0 -~ 4 - 7 -^ 
0.0460 
0. 1590 
_ 0 ._O 441 O- 
0, 1480 
0,0200 

0,0100 

(' . 1 0 9 0 

. O -..0 

0,0600 

0.0 


— il, 2 _ 00 a 

0.9000 i.onoo 

0.7H00 0.9200 

■7-6-1) 0— — ■ — J-. -9-0-0 0- 
2.7400 


0. 1640 
_e_.(V 4 - 74 i. 

0. 1530 
0,0450 
-Ov- 1 ^ 0 - 
0,0430 
0,1300 
- 4 hr 050 - 0 “ 
0.1160 
0.0 1 OO 

O.O 

0 ,093(1 

-O-.-O 

0 ,0800 

0.0 

4.Q-0-0- 
1 . t 0 0 0 
1 .0600 
- 2 -.- 0-^0 0 - 


0,1360 

-0_jCL4-34)- 

0,1240 

0.0400 

0,0370 
0 , 1020 
-O-rO-l-fM)- 

0,0910 

0.0 

-O-s-0- 84)41- 

0.0 

0,0700 

-t)- 4 )- 

0,0600 

0.0 

4L.,-5-0-fi.O- 
1 .2000 
I ,2000 


0,1100 


original page is 

OP POOR QUALM 
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W=604 


Ul 

<J1 



t Numbers betv/een brackets designate the mesh location, as specified in Fig, 3d, 








VD 



ur = 714 ur = 1149 wr = 1319 ur = 1437 

'?] [4, 71 [7, 7] [10, 7] 



ur = 409 
[1, 9] 



'jr =34 3 
[ 1 , 10 ] 



uir — 973 w’r = 1205 wr = 1376 

[4, 9] [7, 91 [10, 9] 



[4, lOi 


[7, 101 


no, 101 


+ Muml-jers betv/een brackets designate the Jnesh location, as specified in Fig. 3d. 


FIG, 3 c, VELOCITY DIAGRAMS (FT/SEC), AT DIFFEKHT MESHES. 
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FIG. 3d. I^IESH coordinates 


Mesh 

Location 
[I, J] 

Radial 
Coordinate 
r, inch 

Axial 
Coordinate 
z , inch 

Mesh. 
Location 
[I, JI 

Radial 
Coordinate 
r, inch 

Aocial 

Coordinate 
z, inch 

[1, 1] 

2.740 

0.0 

[1, 7] 

1.025 

0.4125 

[4, 1] 

2 . 740 

0,153 

[4, 7] 

1.651 

0.685 

[7, 1] 

2,740 

0.305 

[7, 7] 

1.894 

0.792 

I — 1 

O 
* — ! 

2.740 

0.452 

[10, 7] 

2.063 

0.866 

[1, 3] 

2.250 

0.025 

[1, 9] 

0 .587 

0 . 875 

[4, 3] 

2 , 344 

0,232 

[4, 91 

1.397 

0.978 

[7, 3] 

2.419 

0.397 

[7, 9] 

1.731 

1.021 

UO, 3] 

2.481 

0.533 

[10, 9] 

1.976 

1.053 

[1, 5] 

1.549 

0 . 162 

[1, 10] 

0.500 

1.125 

[4, 5] 

1.908 

0 . 451 

[4, 101 

1.278 

1.143 

[7, 5] 

2.091 

0.596 

[7, 10] 

1.665 

1.153 

[10, 5] 

2.215 

0.697 

[10, 13] 

1.950 

1.161 
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Case No, 1 


Case No. 4 


Case No. 



FIG. 5, iMERIDIOfiAL VELOCITY DISTRIBUTION CFT/SEC) 











2600 




Relative velocity (ft/sec) 



FIGURE 8, VELOCITY DISTRIBUTlOll HEAR HUB. 
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FIGURE 14, RELATIVE EXIT FLOW AilGLE (^ 3 ) AND GEOMETRICAL EXIT 
ANGLE (f) )- (TAKEN FROM REF, 7), 


4 8 




FIG, 15a, ILLL'STPv'\TIO:i of relative SL'IFL FCTlOiT, 


4 ? 




FIG, ISb, typical DISTRIBUTIOd FOR THE T,«:iGEilTIAL C0!TO!IEHT 
RELATIVE SWIRL VELOCITY AT ROTOR EXIT. 
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‘■MERIDIONAL PLANE 


FIGURE 16A, NOMEflCLATlJRE OF RELATIVE FLOW IH RADIAL TUPBIHE, 





FIG. 16C: BOUNDARY CONDITIONS FOR THE DISPLACEMENT FLOH 
(RADIAL TURBINE CASE). 
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6g = sin”^(0/t) 



FIC7URE 17A. GEOMETRIC BLADE EXIT AflGLE. 





FIGURE 17E. TYPICAL VELOCITY DIAGRAM AT 
BLADE EXIT (DESIGH POIPT). 
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FIGURE ISA. EFFECT CF SWIRL VELOCITY NEAP SHROUD, 



FIGURE 18B, EFFECT OF SWIRL VELOCITY NEAP HUB, 
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Non-dimensional swirl velocity 





Non-dimensional swirl velocity U 







Non-dimensional swirl velocity 



Radius Ratio, r/r^hroud 


FIG. 23. EFFECT OF CHAflGIIlG NEWER OF BLAPES. 







Relative exit flow anale 



FIG, 25, EFFtCT OF CHAflGIHG THE HASS FLOW RATE O'l THE EXIT AilGLE 


DlSTRIBUTIOn, 
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Relative exit flov; angle y degree 




62 


